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Abstract In this review we confront the current theoretical understanding of particle ac- 
celeration at relativistic outflows with recent observational results on various source classes 
thought to involve such outflows, e.g. gamma-ray bursts, active galactic nuclei, and pulsar 
wind nebulae. We highlight the possible contributions of these sources to ultra-high-energy 
cosmic rays. 

1 Introduction 

High-energy astrophysical phenomena stem from the generation of powerful flows emanat- 
ing from supernova explosions, gamma-ray bursts (GRB), from ejections in the environment 
of black holes or neutron stars that lead to the formation of very strong shocks, at which par- 
ticle acceleration takes place. The new developments in these issues, especially for relativis- 
tic shocks, are based on the interdependence between the shock structure, the generation of 
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supra-thermal particles and the generation of turbulence. It is thought, and numerical simula- 
tions support that view, that the penetration of supra-thermal particles in the shock precursor 
generates magnetic turbulence which in turn provides the scattering process needed for par- 
ticle acceleration through the Fermi process. This successful development, first elaborated 
for supernova remnants (SNR), inspired similar investigations for the termination shock of 
GRBs. However, in ultra-relativistic shocks, difficulties arise with the transverse magnetic 
field that places a limitation to particle penetration upstream and that drags particles in the 
downstream flow and makes shock recrossing difficult. It turns out that only sufficiently 
fast micro-turbulence can make the Fermi process operative, as demonstrated by recent nu- 
merical simulations. Following a review of the main observational results on GRBs, active 
galactic nuclei (AGN), and pulsar wind nebulae (PWN), these points are briefly discussed 
and astrophysical consequences are drawn. We describe the role relativistic shocks inside 
relativistic flows, e.g. the internal shocks of the prompt-emission stage of GRBs, may play 
in the generation of ultra-high-energy cosmic rays (UHECR). Noting that the energy re- 
quired for supplying sub-GZK UHECR is huge compared with the available energy budget, 
we also discuss other sources, such as AGN and young pulsars, that may contribute to the 
flux at ultra-high energies, all the more so if the composition is enriched in heavy nuclei, as 
suggested by recent experimental results. 



2 Gamma Ray Bursts 

The first GRB was observed by one of the Vela satellites monitoring for the Nuclear Test 
Ban Treaty in 1967, but the ( unexpected) astronomic al results were not declassified and 
published for another six years (K lebesadel et al .119731) . For many years the nature of GRBs 
was unknown, since the distance scale was completely unknown. Beginning about 20 years 
ago, the cosmological spatial distribution of GRBs was strongly hinted at due to the very 
isotropic distribution on the sky of GRBs localiz ed by CGRO. The w ealth of detailed infor- 
mation garnered in the last seven years by Swift dGehrels et al.ll2004h has taken the study of 
GRBs to the next level, and indeed the current situation is in some sense more confusing 
than our naive pre-Swift picture (e.g. iGehrels et ap2009h . 

In this subsection we (i) review the basic properties of the two main types of GRBs, 
long and short, (ii) look at long GRBs in more detail, (iii) review the brief history of short 
GRBs and the difficulties entailed in their study, (iv) provide an overview of the acceleration 
processes for GRB jets, and (v) conclude with a recent results on high energy emission 
observed by Fermi. 

2.1 GRB Properties 

GRBs come in two kind s, long and short, where the dividing line between the two is ~2 s 
dKouveliotou et alJll993h . Long GRBs (lGRBs) are thought to be due to the collapse of a 
massive star, while short GRBs (sGRBs) are inferred to be neutron star - neutron star (NS- 
NS) mergers. A further division can be made spectrally according to their hardness ratio 
(i.e., ratio of high to low energies). The redshift range is from about 0.2 to 2 for sGRBs, 
with a mean of about 0.4. For lGRBs the range is between about 0.009 and 8.2, with a mean 
of about 2.3. The typical energy release is ~10 49 — 10 50 erg for sGRBs and ~10 50 — 10 51 erg 
for lGRBs. These ranges are based on observed isotropic-equivalent energies of ~10 51 erg 
for sGRBs and ~10 53 erg for lGRBs, and estimates for jet beaming for each class, 0j ~ 5° 
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for lGRBs and 0j ~ 15° for sGRBs ( iBurrows et al.l2006l : lGrupe et al.l2006h . Beaming angles 
for sGRBs are still highly uncertain. The corresponding beaming factors ft, = 1 — cos 9 S ~ 
df/2 are roughly 1 /300 for lGRBs and 1 /30 for sGRBs. The sGRBs have weaker X-ray 
afterglows, a mean value of ~7xl0~ 10 erg cm~ 2 s _1 versus ~3xl0~ 9 erg cm -2 s _1 for 
lGRBs. Figure Q] shows spectra for several representative GRBs, and two other high-energy 
sources, the Crab nebula and Cyg X-l. 
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Fig. 1 Representative broad-band vFy spec tra jGehrels et al.l2009l) of a 1GRB ( 910503i llKaneko et all2 007) 
and a sGRB (980425i iKaneko et alj|2003> along with the Crab pulsar nebula iKuip er et al. 200l|) and Cyg 
X-l jMcConnell et alj2002l) . 



2.2 LGRBs 

The BeppoSAX m ission made the critical discovery of X-ray afterglows of long bursts 
dCosta et al] 1997 ). With the a ccompanyi ng discoveries by ground-based telescopes of op- 
tical i van Paradiis et al.lll997h and radio dFrail et al.lll997h afterglows, lGRBs were found 
to emanate from star forming regions in host galaxies at typical distance of z — 1 — 2. Bep- 
poSAX and the following HETE-2 mission also found evidence of associations of GRBs with 
Type Ic SNe. This supported the growing evidence that lGRBs are caused by "collapsars" 
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where the central core of a massive star collapses to a black hole jMacFadven & WooslevI 
Il999h . 

LGRBs are incredibly bright. A typical galaxy at a redshift of only z = 3 is fainter than 
m ~ 27. Multiwavelength observations of the current record holder, GRB 090432 (at z — 8), 
are providing information about the universe at a time when it was only abo ut 4% of its 
current age, and shed light o n the process of reionization in the early universe dTanvir et al.l 
120091 : ISalvaterra et alj|2009r) . The highest redshift GRBs are seen to have high luminosity, 
resulting in fluxes well above the detection threshold. Such bursts are also strong at other 
wavelengths. Table [T| presents optical data for the highest redshift GRBs observed to date, 
where the look-back time /LB(Gyr) is given in column 2. 



Table 1 High z GRBs. 
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2.3 SGRBs 

At Swift's launch, the greatest mystery of GRB astronomy was the nature of short-duration, 
hard-spectrum bursts. Although more than 50 lGRBs had afterglow detections, no afterglow 
had been found for any sGRB. Swift provided the first sGRB X-ray afterglow localization 
with GRB 050509B and HETE-2 enabled the first optical afterglow detection with GRB 
050709. These two bursts, plus Swift observations of GRB 050824, led to a breakthrough i n 
our understanding jGehrels et al .1120051 |B loom et al .1120061: iFox et alj|2005l : Ivillasenor et all 



l2005l : lHiorth et al.ll2005l : lBarthelmv et al.l2005klBerger et al.l2005l) of sGRBs. BAT has now 



detected 60 sGRBs, most of which have XRT detections, and about one third of which have 
host identifications or redshifts (an additional two have been detected by HETE-2, one by 
INTEGRAL, and two by FermifLAI). We now have ~50 sGRB localizations. 

In stark contrast to long bursts, the evidence to date on short bursts is that they can 
originate from regions with low star formation rate. GRB 050509B and 050724 were from 
elliptical galaxies with low current star formation rates while GRB 050709 was from a 
region of a star forming galaxy with no nebulosity or evidence of recent star formation ac- 
tivity in that location. Recent HST observations of locations of sGRBs in their hosts reveal 
that sGRBs trace the light distribu tion of their hosts while lGRBs are concentrated in the 
brightest regions dFong et al.ll2010h . SGRBs are also different from lGRBs in that accom- 
panying supernova e are not detected for nearby events dBloom et al ] l2006l :l Fox et al ] l2005l : 
iHiorth et al.ll2005h . Taken together, these results support the interpretation that short bursts 
are associated with an old stellar population, and may arise from mergers of compact bina- 
ries [i.e., double neutron star or neutron star - black hole (NS-BH) binaries]. 
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2.4 GeV Emission 

Fermi was launched into low-Earth orbit in June 2008 and has two primary high-energy 
detectors: the Large Area Telescope (LAT) which operates between 20 MeV and > 300 GeV, 
and the Gamma-Ray Burst Monitor (GBM) which operates between 8 keV and 40 MeV. So 
far the LAT has detected 24 GRBs; two were sGRBs, and nine showed extended emission. 
The emission from GRB 090902 B included a 34 G eV photon. One of the most luminous 
to date has been GRB 080916C i lAbdo et al.|[2009ah at a redshift of 4.35. It had extended 
emission (18 min) and exhibited a lag in LAT energies with respect to GBM. GRB 090510 
is unique in being the only short burst with a known redshift (0.903) showing GeV emission. 
The lack of detectable time delay between specific peaks in the light cu rves of GRB 09051 
at different energies lead s to strong constraints on Lore ntz invariance ( lAbdo et al.ll2009bh . 
Recent theoretical work dKumar & Barniol Duranll2010l) on the Fermi /LAT detected GRBs 
suggests that these may represent unusually powerful explosions with Lorentz factors >10 3 
in which the entire progenitor is obliterated. The simplest model, namely an external shock 
with synchrotron emission, can be used to take the early values (at ~10— 10 2 s) of the 
observed high-energy emission and successfully predict the much later values of the optical 
and X-ray afterglow (at ~10 5 - 10 6 s). 



2.5 GRB Summary 

Recent progress in GRB research has been strongly motivated by observational discoveries. 
To date, Swift has detected about 600 GRBs, FermifLPS 24. High redshift GRBs are illumi- 
nating the properties of the high— z universe and probing into the era of re-ionization. Swift 
finds sGRBs in different environments than lGRBs; also sGRBs are not accompanied by su- 
pernovae. The accumulating evidence provides support for the NS-NS merger model. Many 
GRBs have delayed onset of GeV emission, and more have extended high energy emission. 
Interesting constraints on the Lorentz factors associated with outflow, and Lorentz invari- 
ance violation, also come from the synergism between GeV and lower energy observations. 
It is not currently know for certain whether GRB jets are made primarily of baryons or 
Poynting flux, but momentum for the latter idea appears to be gaining strength. 



2.6 Jet Launching Processes 

There are currently two primary lines of thought regarding the creation and propagation of 
jets in GRBs. 

The baryonic jet model, whose roots can be traced back to t he elegant analytical solu- 
tions of a relativistic blast wave by iBlandford & McKeel ( fl976h . posits that a jet contain- 
ing about a Jupiter's mass worth of gas, ~ O.OOIMq, is som ehow launched nea r the BH 
created by the collapsar with a Lorentz factor F ~ 10 — 20 dzhang etal.ll2003h . The jet 
propagates through the dense stellar envelope of the progenitor star where it is focussed 
and compressionally heated. After breaking free of the stellar surface, the thermal energy 
of the compressed jet is translated into bulk kinetic motion, further accelerating the jet to 
a Lorentz factor of several hundred. Subsequent deceleration by the circumstellar medium 
of the jet, which is then idealized as being the fragm ent of a relativi stic shell so that the 
Blandford & McKee formalism can be brought to bear dSari et alj 19981) , can then be used to 
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infer jet beaming angles from putative achromatic "breaks" in the GRB decay light curves 
dFrail et alfcOOlh . 

The Poyntin g flux jet model, which has been gaining momentum in recent years, traces 
its roots back to Blandford & Znaieklll 19771) who considered the electromagnetic (EM) ex- 
traction of energy from within the ergosphere of a Kerr BH. In the last few years workers 
have developed sophisticated numerical codes that calculate the 3D evolution of gas and 
EM fields from the inner edge of accretion disks onto spinning BHs, taking into account 
both general r elativity and magnetohydrodynamics (IMcKinnevl2005l : [McKinnev & Naravanl 
I2007all2007bh . The accretion disk inner edge provides a natural collimating surface. These 
workers find a baryonic zone-of-exclusion within the jet, which effectively suppresses any 
baryonic component. For numerical stability the numerical codes need to have some mass 
within each grid point, therefore a small trace amount of matter is constantly added within 
grid points where the EM fields try to exclude it. These are called "floor" models. Previous 
studies that attempted to place constraints on the jet Lorentz factors, which were based on 
the baryonic jet assumption, have been called into question. The ramifications on the jet 
dynamics of having a predominantly Poynting flux jet have not been developed yet in any 
detail. 



3 Nonthermal Particle Acceleration in GRBs: challenges and perspectives 

Physical mechanisms of prompt emission in GRBs are still to be established. There are fun- 
damental questions of the powerful relativistic outflow composition and matter vs. magnetic 
field dominance to be addressed with both observations and advanced models. The relativis- 
tic outflows may be different for two types of GRB progenitors under consideration. Long 
and soft GRBs are most likely connected to supernovae, while the short and hard GRBs are 
possibly related to compact binary mergers. 

Magnetized outflows from GRB engines of different nature have been studied for some 
decades (see. e.g. | Usovll992l : lThompsonll994l : lMeszaros and Reesll997l : lLvutikov and Blandford! 



l2004l ; lzhang and Yanll2011l ). The outflows converting the rotation power of a compact col- 
lapsar into a broad band of radiation are likely models of GRBs and AGNs. The principal 
question here is our understanding of the microscopic mechanisms of the conversion of 
magnetic field energy into non-thermal particles and the observed emission - with or with- 
out shock formation. Anisotropic striped wind with alternating magnetic polarity is consid- 
ered a favorable configur ation to convert the magnetic energy into the observed radiation of 
pulsar wind nebulae (e.g. iKirk et al ] l2009l :l Sironi & Spitkovskv 201 1). An alternating mag- 
netic field configuration may occur in relativistic jets of GRBs. Magnetic field reconnection 
demonstrated to produce electric fields that accelerate particles in the Earth magnetosp here, 



solar flares, and some laboratory plasma config urations (see e.g. Yamada et al. 2010h 



also proposed as a plausible GRB model by (e.g.jSpruit et al.l200ll:lLvutikov and Blandford! 
|2004| ; iGiannios and SpraitH2005l ; Izhang and Yanlboill i lMcKinnev and Uzdenskvll2012l) . 

The most elaborated model of the origin of the GRB prompt emission is, by now, 
the relativistic dissipative fireball model. There are a number of alternative particle ac- 
cele ration and radiation processes within t he relativistic dissipative fireball paradig m (see 



e.g. iRees and Meszaroslll994 |Piranll2004 lMeszarosll2006l : iMimica and Aloy||2012r) . Syn 



chrotron and inverse Compton radiation in the optically thin regions of the relativistic fireball 
can be associated with the non-thermal electron/positron accelerated either in the internal 
dissipation processes due to shocks or to the flow magnetic field reconnections. Another po- 
tential component is the photospheric emission (thermal or non-thermal) that is rather rarely 
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Fig. 2 Simulated energy dissipation efficiency rj denned by Eq|5]as a function of the shell magnetization 
parameter o~ for the inelastic collision of two initially cold magnetized shells of equal masses, but of different 
Lorentz factors (F = 500 and F2 = 1,000). The key parameter is the final magnetization of the hot merged 
shells Of was chosen to be 0.1 (dot-dashed curves), 1.0 (solid lines) and 10.0 (dashed lines). The two curves 
of the same line style in the Figure differ by the assumed initial magnetization parameters. The top curve for 
each line style corresponds to the case of equal initial magnetization C7[ = o~2, while the lower curves of each 
type correspond to the fixed o~2 = 0.1. The adiabatic index 7 = 4/3 was fixed for the matter in the hot merged 
shell. 

identified in the observed GRB emission. The prompt emission light curves and spectra are 
generally in agreement with the internal shock models. However, some potential problems 
of the scenario are its efficiency and the lack of a bright photospheric component observed in 
a few GRBs. We shall discuss now some general features of the internal dissipation models 
with emphasis on particle acceleration processes. 

3.1 Energetic Efficiency of Internal Dissipation Models: Shell Collisions in Jets 

An important issue of the internal dissipation scenario where the energy carried out by mul- 
tiple colliding shells of different magnetization is the efficiency of the conversion of the out- 
flow power into the observed radiation. The variability of the central engine (of a timescale 
t vai ) driving the relativistic outflows of the mean Lorentz factor F, either matter or Poynt- 
ing dominated, can be modelled as a collision of energy-containing shells. T he model is 
considered to explain the main features of the GRB prompt emission (see e.g. |Piranll2004l ; 



g 




Fig. 3 Simulated dimensionless mean enthalpy that is characterizing the mean Lorentz factor of the random- 
ized particles in the rest frame of a hot merged shell. The curves are simulated for the same parameter sets as 
it is indicated in Figure|2] 



lMeszarosll200fjh , as a vital alternative to the photospheric models of GRBs. The dissipation 
region typically exists at the radii about rdj ss ~ cf val -F 2 . In the case of matter dominated 
jets (of low magnetization) the inner dissipation occurs in the internal shocks while in the 
electromagnetically (Poynting flux) dominated jets the magnetic field reconnection effects 
are most likely crucial though shocks may also occur. The microphysics of the dissipa- 
tion in relativistic s hocks as well as modeling of the magnetic field r econnections are un- 
der intense studies I Yamada et al. 20 id : iBvkov and Treumannl l201lL ISironi & Spifkovskvl 
1201 ll : iMcKinnev and Uzdenskvl 1 20121) 7 Realistic models of a jet that would simulate the 
global RMHD dynamics and simultaneously resolve the dissipative microphysical plasma 
processes at much smaller spatial scales are not feasible at the moment. However, simple 
multiple shell models of the internal dissipation that just parameterize the magnetic field re- 
connection effects are still rather useful and it is instructive to discuss some of these models. 
A similar approach can be applied to other relativistic outflows, like those of the AGN jets 
(see for a discussion <j4j and of the pulsar wind nebulae. 

To illustrate the effect of the outflow magnetization o~ = B 2 /Anrpc 1 on the energy con- 
version into the observed radiation it is instructive to use a simple two shell model (see e.g . 
iPanaitescu et al ]|!999l ; lKumaj|l999l ; lDaigne and Mochkovitchll 19981 ; Izhang and Yanl,201lh . 
The shell collision may result in dissipation of the magnetic energy due to reconnection 
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Fig. 4 Simulated Lorentz factor if of a hot merged shell in the observer frame. The curves are simulated for 
the same parameter sets as it is indicated in Figure [2] 



and turbulence cascade. Consider two shells of masses, Lorentz factors and magnetization 
parameters [((mi, 1], Of) and (m 2 ,r 2 ), 2 )], respectively, colliding inelastically with the for- 
mation of a merged shell of (mf,If, Of). The internal energy 8' released in the rest frame of 
the merged shell (of the Lorentz factor if) is assumed to be in the form of either thermal 
or non-thermal accelerated particles and radiation with nearly isotropic distribution in the 
rest frame. Then, in the observer frame the released energy is AE = if 8'. This leads to con- 
version of some amount of magnetic energy into internal energy of the fluid, and then to 
radiation. Let us envisage a picture where the two shells merge with a lower magnetization 
parameter Of by the end of such an inelastic collision. Energy conservation and momentum 
conservation can be presented as 

r l W l m l +r 2 W 2 m2=r f Wf(m l +m 2 + 8')-^-8', «ff(JT) = 1 + ^fc-U (1) 

and 

r,ftmi ( 1 + 01 ) +r 2J S 2 m 2 ( 1 + o 2 ) = r f /3 f (mi + m 2 + 8') ( 1 + o f ) . (2) 
Since the merged shell can not be considered as a cold one anymore 

(3) 
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where is the proper kinetic energy density (i.e. the internal energy density with the rest 
energy density subtracted), V' is the shell volume, and P' is the pressure in the rest frame of 
the merged shell. Note that to calculate the energy released in the rest frame of the merged 
shell 8' one should keep all of the terms of the order of F ~ 2 . The expression for f{(i7) ls 
exact for the case of the transverse magnetic field in the rest frame of a shell, and it accounts 
for the energy of the induced electric fields in the observer frame. 

In general, there are a few distinct components that contribute to the pressure and the 
proper energy density (thermal and nonthermal baryons, leptons, and photons). If for sim- 
plicity we describe these as a single fluid with adiabatic index 7, then 

and one can resolve the energy-momentum conservation equations, assuming the simple 
equation of state Eq|4]to calculate the Lorentz factor of the merged shell fj and the energy 
dissipation efficiency r/ of the inelastic collision of the two cold magnetized shells 

= rfS ' (5) 

il fimi +r 2 x I , 2in2 

In Figure [2] we show the simulated energy dissipation efficiency r\ as a function of the 
shell magnetization parameter o for the inelastic collision of two initially cold magnetized 
shells of equal masses, but of different Lorentz factors. We consider the GRB jet as a generic 
case and therefore choose T\ = 500 and i~2 = 1 , 000. The energy dissipation efficiency 77 , the 
mean enthalpy of the hot matter in the merged shell, that characterizes the mean Lorentz 
factor of the randomized particles in the rest frame of the hot merged shell (shown in Figure 
[3), and the Lorentz factor if of the hot merged shell in the observer frame, shown in Figure 
[4] are derived from the conservation laws and the equation of state. The adiabatic index is 
fixed to 7 = 4/3, though in more accurate numerical simulations it depends on the the the 
Lorentz factor If. The key parameter here is the final magnetization of the hot merged shell 
Of that is determined by the currently poorly known rate of magnetic field reconnection 
in the merging shells. The case of fast field reconnection (and, therefore, efficient magnetic 
field dissipation) is illustrated by the dot-dashed curves in Figure [2] where Of = 0. 1 . The cases 
of lower magnetic field dissipation are illustrated by the final magnetization parameters Of 
= 1.0 (solid lines) and Of = 10.0 (dashed lines). 

The dependence of the merged shell parameters on the initial magnetization of the col- 
liding shells is presented by two curves of the same style that differ by the initial magne- 
tization parameters. The top curve for each type of line corresponds to colliding shells of 
equal initial magnetization 01 = 02. The lower curves of each type correspond to the the 
case when the fast shell of the Lorentz factor Ij = 1 , 000 has low initial magnetization 02 = 
0.1. 

It is clearly seen in Figure [2] that the dissipation efficiency is higher in the case of the 
initially highly magnetized shells with fast magnetic field reconnection resulted in low Of 
of the merged shell (dot-dashed curves). Again, the microscopical model of the reconnec- 
tion rate in the complex flow is still to be done to estima te the crucial value of the final 
magnetization parameter Of of the merged shell (see e.g. lYamada et al. ||2Ql3). To inves- 
tigate the problem of shock formation in the internal dissipation scenario, a microscopic 
modeling of the collisionless shock formation and its structure in a h ighly magnetized rela- 
tivistic outflow is needed, and tha t is a truly challenging task (see e.g. lBvkov and Treumannl 
1201 ll : ISironi & Spitkovskvl 1201 ll) . The problem of Fermi acceleration in transverse rela- 
tivistic shocks of different magnetization that are important to describe the external shocks 
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in the jets of GRBs and AGNs will be addressed in ij6] here we concentrate on the internal 
dissipation models. 

As it is seen in Figure [4] the high Lorentz factors if > 1 , 000 of a hot merged shell in 
the observer frame can be achieved even for the incomplete magnetic field dissipation in the 
merged shell (solid and dashed curves) and, therefore, it can further catch up other slowly 
moving shells and merge with them providing a chain dissipation process (Bykov & Osipov 
2012, in press). 

The mean Lorentz factors in Figure |3]derived from the conservation laws do not preclude 
a presence of a non-thermal (e.g., a piece-wise power-law) particle distribution, where some 
minor fraction of particles can reach energies that are by some orders of magnitude larger 
then the derived "thermodynamic" mean Lorentz factor. Now we turn to discuss in brief 
possible particle acceleration processes. 



3.2 Particle Acceleration in the Internal Dissipation Models 

Relativistic turbulence produced by the internal dissipation in shocks and magnetic re- 
connection in relativistic jets result in acceleration processes occuring on both gyro time 
scale and on longer comoving hydrodynamical time scales of the order of l/c. The elec- 
tric fields induced by turbulent motions of plasmas carrying magnetic fields on different 
scales lead to statis tical energy gains of the superthermal charged particles and their wide- 
band radiation (see iBvkov and Meszaros 19961: iMizuno et al]|201lLlNishikawa et al. 20 id : 



ISironi & Spitkovsky||201ll ; lzhang and YarfeOl lllBvkov et alfcouJlVIurase et alj|2012l) 



For nonrelativistic MHD turbulence the particle energy gain over a turbulent correla- 
tion length (or correlation time) is small, because the induced electric field is smaller then 
the entrained magnetic field. However, the distinctive feature of statistical acceleration in 
the relativistic MHD turbulence and shocks on larger scales expected in the flow-colliding 
regions, is the possibility of a substantial particle energy gain over one correlation scale, be- 
cause the induced electric fields are no longer small. In this case the standard Fokker-Planck 
approach cannot be used. Instead, it is possible to calculate the energy spectra of nonther- 
mal particles within a specia l integro-differential equation which is a generalization of the 
Fokker-Planck approach (see lBvkov and Toptygin|[l993l : lBvkov and Meszaroslll996h . 

Charged particles interact with a wide spectrum of RMHD fields and an internal shock 
ensemble produced by the colliding shells. In the comoving frame, it is assumed that the fluc- 
tuations on all scales up to ~ A (including the internal shock ensemble) are nearly isotropic 
(in the latter case, it is enough if they are forward-backward symmetric). The small mean 
free path A of the superthermal particles leads to their isotropy in the frame of the local bulk 
velocity fluctuations. The assumed statistical isotropy of the bulk velocity fluctuations in the 
comoving frame of the wind results then in a nearly isotropical particle distribution, after 
averaging over the ensemble of strong fluctuations on scales ~ I. 

To calculate the spectrum of nonthermal leptons accelerated by an ensemble of internal 
shocks and large-scale plasma motions in the flow-colliding region (averaged over the sta- 
tistical ensemble of large-scale field fluctuations) we use a kinetic equation for the nearly- 
isotropic distribution funct ion N = y^F, which takes into account the non-Fokke r-Planck 
behavior of the system (see lBvkov and Toptygin|[l"993l ; lBvkov and Meszarosll 19961) . 

= &(«)+n,d&Z> 1 (§ - ^AFbSut) 
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Here = ln(y/y), y is the Lorentz factor of the injected particles, °= £c«/ 2 is the rate 

of nonthermal particle injection, n is the lepton number density in the local flow comoving 
frame. The kernels of the integral equation Eq.(|6} determining the spatial and momentum 
diffusion are expressed through correlation functions which describe the statistical proper- 
ties of the large scale MHD turbulence and the shock ensemble. Following the renormal- 
ization approach, the Fourier transforms of the kernels Df (s) and D 2 is) are solutions of a 
transcendental algebraic system of equations of the form D^ 2 — ^i.2(Df ,-Df ,s). Here s is a 
variable which is the Fourier conjugate of . Equation $6$ and the renormalization equations 
are valid only for particles with sufficiently small mean free paths A(y) <C A. 

It is important to note that the solution of equation © has a universal behavior, only 
weakly dependent on the complicated details of the turbulent system. The stationary solution 
to Eq.© with a monoenergetic injection rate Qi has an asymptotical behavior of a power-law 
form, N oc g,y~ a , where a = -0.5 + [2.25 + 6D1 (0)L>2 1 (0)] 5 , and thus one may take 6 ~ 
(l/A) 2 . For conditions typical of a developed RMHD turbulence, the ratio of the rate of the 
scatterings to the acceleration rate is D\ (0)0^ 1 (0) < 1, and for 8 < 1 one obtains a ~ 1. This 
hard y~ l spectral behavior arises because the acceleration time X a ~ l/c ~ aA/c is much 
shorter than the escape time at the relevant energies, % esc ~ A 2 /k ~ A 2 /(lc) ~ 4 /(cue). The 
power needed to produce such a spectrum of nonthermal particles increases °c y max , so it is 
important to understand its temporal evolution. 

In the test particle limit, where the backreaction of the accelerated leptons onto the 
energy-containing bulk motions is small, we have N(y,t) £«y~' for y< y+(f), where 
Y*(t) = yiexp(t/i%)tmd 

T h a <*l/c~a{A/c), (7) 

is the typical hydrodynamical acceleration timescale (see, e.g. iBvkov and Toptygin|[l993h . 
with Yi ~ few, a < 1 , and the comoving width of the region energized by shocks equal to A . 
From the energy balance equation, when the value y*(f) ~ y,- m p /m e e £ is reached, the 
growth must saturate, and the resulting spectrum consists of two branches. One of these is 
the hard spectrum N(y) ~ £ny ~' , for y < y,, where 

y*~ ym p /m e £ . (8) 

For the typical scales of the considered problem y ~ 1 and C, ~ 10 -3 , so y + ~ 10 5 (but it 
could be even larger since e ~ 1 for large scale plasma motions). In the regime of very effi- 
cient particle acceleration where the backreaction of the accelerated leptons on the energy- 
containing bulk motions is substantial, one should use a non-lin ear approach . In that case the 
partic le distribution experience fast temporal evolution (see e.g. lBvkovl200ll ; lFerrand and Marcowithl 
l2010h . 

The multiwavelength observations discussed above can be used to constrain the charac- 
teristics of relativistic turbulence and particle spectra evolution. Measurements of polariza- 
tion of prompt gamma-ray emission of GRBs would be very informative. Polarization mea- 
surements provide an additional deep insight into the nature of acceleration/radiation pro- 
cesses. By now only a few instruments are available for this kind of measurement at gamma- 
ray energies. Gamma-ray bursts (GRBs) are very promising candidates for polarimetry due 
to their large flux over limited time intervals, maximizing the available signal-to-noise ratio. 
To date, a few polarization measurements have been reported, all claiming a high degree 
of polarization in the prompt emission of GRBs, but with rather low statistical evidence 
vet. iGotz etaD J2009h used the IBIS telescope aboard the INTEGRAL space laboratory to 
measure the polarization of the prompt gamma-ray emission of the long and bright GRB 
041219A in the 200-800 keV energy band. They found a variable degree of polarization 
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ranging from less than 4% over the first peak to 43%-25% for the whole second peak. Time- 
resolved analysis of both peaks indicated a high degree of polarization, and the small av- 
erage polarization in the first peak can be ex plained by t he rap id variations observed in the 
polarization angle and degree. The results bv lGotz et al. I (l2009n are consistent with different 
models for the prompt emission of GRBs at these energies, but they favor synchrotron radi- 
ation from a relativistic outflow with a magnetic field, which is c oherent on an angular size 
comparable with the angular size of the emitting region. Recently. I Yonetoku et all 1 1201 lh re- 
ported the polarization measurement of the prompt gamma-ray emission of GRB 100826A 
with the Gamma-Ray Burst Polarimeter (GAP) aboard the small solar-sail demonstrator 
IKAROS. They detected the change of polarization angle during the prompt emission, and 
the average polarization degree of 27± 11%. Polarization measurements are a powerful tool 
to constrain the GRB radiation mechanisms and the magnetic field structure in the photon 
emitting regions. 

4 Active Galactic Nuclei 

4. 1 Studies of Particle Acceleration in AGN Jets 

The relativistic outflows of Active Galactic Nuclei (AGNs) are well-studied particle acceler- 
ators. The AGN jet phenomenon spans many orders of magnitudes: jets originate on sub-pc 
scales close to supermassive black holes with typical masses between a few million and a 
few billion solar masses and can propagate over many hundred kpc to feed giant hot-spot 
complexes and radio lobes. One of the most remarkable properties of jets is that they dis- 
sipate little energy while transporting vast amounts of energy and momentum over large 
distances. However, jets are not dissipation-free neither one sub-pc scales nor on kpc-scales: 
they do accelerate electrons and maybe also protons to high energies. The high-energy par- 
ticles interact with magnetic fields, photon fields, and/or ambient matter and make AGNs 
some of the brightest extragalactic sources of continuum emission across the electromag- 
netic spectrum. There are many reasons to study particle acceleration in AGNs: we would 
like to explain the observed electromagnetic radiation in order to constrain the composi- 
tion and structure of jets. The composition and structure constrain the processes of AGN 
accretion and jet formation, acceleration and collimation. The same studies can enhance our 
understanding of AGN feedback, i. e. how AGNs interact with their hosts and decisively 
impact their evolution. Studies of AGNs also allow us to study the particle acceleration 
mechanisms. A good (microscopic) understanding of the dissipation processes in jets is in- 
strumental for addressing the larger questions mentioned before. 

In the following we discuss two particle acceleration sites: the blazar zone less than a 
pc away from the supermassive black hole, and jet particle acceleration by kpc-scale jets. 
The processes at work in these two well-separated regions have recently received a lot of 
attention mainly because of spectacular jet images from the VLB A, HST, and Chandra, and 
spectacular multiwavelength observations with the RXTE, Suzaku, and Swift, X-ray, and 
Fermi, H.E.S.S., MAGIC, and VERITAS /-ray telescopes. 

4.2 Studies of Particle Acceleration Processes in the Blazar Zone 

The cores of the "blazar class" of AGNs are bright sources of spatially unresolved continuum 
emission. In the case of BL Lac objects, emission and absorption lines are absent or weak. In 
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the case of their more powerful siblings, Flat Spectrum Radio Quasars (FSRQs), the spectra 
exhibit emission and absorption lines, which afford additional diagnostics concerning the 
mass of the central black hole, the accretion power, and the redshift of the source. Blazars 
are sources with jets closely aligned with the line of sight. The relativistic propagation of 
the jet plasma and the close alignment of the jet with the line of sight lead to relativistic 
boosting of the emission. Blazar Spectral Energy Distributions (SEDs) show evidence for 
two broad peaks presumably coming from synchrotron and inverse Compton emission from 
the same electron population. O bservations of the BL Lac Mrk 501 revealed emission up to 
16 TeV jAharonianetal]|l999t) . giving direct proof that AGN jets accelerate particles to 
3> 1 TeV energies. AG Ns may accelerate partic les to much hi gher energies, maybe e ven to 
ultra-high energies (e.g. iKachelriefi et al.ll2010h . see however lLemoine & Waxmanl j2009) 
for a detailed discussion of acceleration power in these sources. 

Two scenarios are commonly invoked to explain the observed emission: (i) the jet is 
initi ally Poynting flux dom inated and accelerates particles in magnetic reconnection events 
(e.g. lGiannios et al1l201(]| . and references therein); (ii) the jet is either particle energy dom- 
inated right from the start, or, it is initially Poynting flux dominated and some unknown 
mechanism converts the energy from Poynting flux into bulk motion energy, and the particle 
dominated jets subsequently accelerate particles at shocks. 

Blazar studies of particle acceleration benefit from the fact that the broad-band SEDs oft 
he sources evolve on short time scales: X-ray and gamma-ray flares with 1 min durations 
have been reported. It thus is possible to track the temporal evolution of the particle en- 
ergy spectra. Modeling of sequences of broad-band SEDs has sh own that the jets are highl y 
relativistic with bulk Lorentz factors of ~50 or even higher (e.g. lKrawczvnski et al.ll200~lh . 
The simplest types of models synchrotron self-Compton models - indicate that the particle 
ener gy density dominates strongly over the magne tic field energy density in the blazar zone 
(e.g. lKrawczvnski et alj|200l : I Acciari et al.ll201lh . These results clearly favor shock accel- 
eration over magnetic reconnection. The data and theoretical studies (particle-in-cell (PIC) 
simulations and instability analyses) favor mildly relativistic shocks internal to the jets over 
highly relativistic external shocks: the former have faster downstream plasmas commensu- 
rate with the high inferred bulk Lorentz factors of the emitting plasma, and are less suscepti- 
ble to the strong suppression of the formation of upstream scattering c enters by even a weak 
plasma magnetization than their highly relativistic counterparts (e.g. ILemoine & Peiletieri 
l20ld : ISironi & Spitkovskvll2009h . 

X-ray and very high-energy y-ray observatories have recently succeeded to sample the 
temporal evolution of the X-ray and -ray energy spectra with delicate accuracy. Although 
observers organized a large number of observation campaigns with the objective to find 
the flux vs. spec tral index evolution patterns predicted by simple accelera tion theory (e.g. 
Kirk et al. 19981), the ob servations revealed rather erratic evolutions (e.g. iTakahashi et alj 



2000l : Garson et af1l2010l and references therein). Unfortunately, observations have not yet 



succeeded to determine unambiguously where exactly the blazar emission originates. For the 
radio galaxy M87 possibly a misaligned blazar the observation of temporally coincident 
radio, X-ray, and y-ray flare indicat es that the emission comes from <100 Schwarzschild 
radii of the supermassive black hole dAcciari et al1l2009b . 



4.3 Studies of Particle Acceleration Processes of kpc-scale Radio, Optical, and X-Ray Jets 



The VLBA, HST, and Chandra telescopes are delivering images of the kpc-jets of an ever 
increasing number of radio galaxies. The images in the radio, optical, and X-ray bands can 
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be used to infer complementary information about the emitting particle populations. Radio 
and optical polarimetry give additional clues about the orientation of the magnetic field in 
the bright jet regions. 

A recent somewhat surprising result was the detection of a large nu mber of X-ray 
bright kpc-jets with Chandra (see the review bv lHarris & Krawczv nski 2006, and references 
therein). In the case of powerful Fanaroff-Riley (FR) type II sources (like PKS 0637-752) 
the combined radio, optical and X-ray energy spectra showed unambiguous evidence for 
two distinct emission components. Presently two models are commonly invoked to explain 
the second X-ray brig ht component: (i) inverse Compton scatte ring of Cosmic Microwave 
Background photons jTavecchio et al]|2000l : [Celotti et al.ll200ll) . (ii) a second synchrotron 
component dHarris & Krawczvnskill2002l) . The first model requires jet plasma moving with 
large bulk Lorentz factors (F ~10) at kpc-distances from the central engine. The relativistic 
motion boosts the CMB photon energy density in the reference frame of the emitting plasma 
by a factor of F 2 , and the mean photon energy by a factor of F. These two effects make it 
possible to explain the observed X-ray emission with electrons with modest electron Lorentz 
factors on the order of y = 100. The second model requires much higher Lorentz factors: 
assuming equipartition magnetic fields, electrons with Lorentz factors y ~ 10 7 are needed to 
explain the X-ray emission. Both models have severe shortcomings. For example, in the first 
model the long cooling times of the low- y electrons are at odds with the well-defined knotty 
structure of some X-ray bright jets. A weakness of the second model is that there is not yet 
a good explanation for the existence of two distinct electron populations with very differ- 
ent spectral properties. For lower-power FR-I-type sources like M 87, Cen A, or 3C 31, the 
combined radio, optical and X-ray energy spectra are consistent with a synchrotron-origin 
of the X-ray emission from a single population of electrons. For an equipartition magnetic 
fields of B ~100 fiG, the X-ray emitting electrons have Lorentz factors y ~ 10 7 , and ra- 
diative synchrotron cooling times of a few years. The X-ray emitting electrons thus do not 
move far from their acceleration sites before they loose their energy, and the X-ray bright- 
ness profiles tracks the acceleration of the high-energy electrons. The diffuse appearance of 
some jets implies quasi-continuous acceleration. 

A few radio galaxy jets have been studied with the Hubble Space Telescope giving not 
only high-resolution images of the optical brightness but also of the optical polarization. 
The magnetic field probed by the optically emitting electrons seems to be aligned parallel 
to the jet flow for the most part. However, upstream of the brightness maxima, the field are 
perpendicular to the jet flow. The radio and optical polarization behavior dif fers, indicating 
that th e emission at different wavelengths samples different regions of the jet l lPerlman et al.l 
Il999h . 



5 Acceleration to ultra-high-energies 

5.1 Some properties of UHECRs 

One of the 11 fundamental s cience questions fo r the 21st century listed in the final report 
of the 2002 Decadal Review ^Turner et al1l2002l) is the nature of cosmic rays. The detection 
of cosmic rays at ultra-high energies (UHECR) dates back to the early Sixties, but only 
during the last 20 years detectors of sufficiently large si ze have become opera tional, that 
the origin of UHECRs can be addressed (for a review see Kotera & O linto 2011). UHECRs 
are indirectly detected by observing the airshowers they trigger in the atmosphere. One 
can look for fluorescence emission or other radiation produced high in the atmosphere, or, 
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alternatively, one registers the passage of secondary particle in charged-particle detectors 
on the ground. Th e AUGER observatory (I Abraham et al.lE004h and the Telescope Array 
l lKawai et alkOQSl) combine both techniques. 

The main observables used to infer the properties of UHECRs are the anisotropy, the 
composition, and the spectrum. The composition is difficult to determine, because the ap- 
pearance of giant airshowers can only be modeled with particle-physics event generators 
that involve extrapolations of behaviour observed in accelerator experiments, for the CoM 
energy of an arbitrary nucleus of 10 EeV energy with a nitrogen nucleus at rest is far higher 
than that achievable with even the largest man-made accelerator, the LHC at CERN. Con- 
siderable systematic uncertainty thus overshadows attempts to study the composition. It ap- 
pears that between 1 PeV and 0. 1 EeV we observe a trend from a pre dominance of light 
particles to heavy nuc lei. Around 1 EeV, the composition is light again jAbbasi etal 
lAbraham et al ■l2010bl) . Above 1 EeV, Auger observes a transition to heavier particles that is 
not see with other experiments at this time, possibly on account of statistics. 

The anisotropy is low around 1 EeV, where upp er limits near 1% have been published 
for the sidereal dipole anisotropy dAbreu et alfeOllh . At higher energies above 57 EeV, for 
which little deflection would be expected, if the primary particles were protons, a correlation 
is observed between the arrival direction of particles and certain types of nea rby AGN which 
in the end are proxies o f the matter distribution within ~75 Mpc from us llAbraham et al.l 
2008l: lAbreu et al.l2010h . The distance limitation is expected because nuclei at these energies 
undergo photodisintegration and photomeson production that provide losses on correspond- 
ing time scales, leading to the so-called GZK cut off. 

Whereas cosmic rays approximately obey a power-law spectrum with index s ~ 2.7 
(dN/dE <x E~ s ) below the so-called knee in the spectrum at 3 PeV, the spectrum of UHE CRs 
is soft between a few PeV and 3 EeV with a power-law index s ~ 3 dApel et al.ll20"09h . At 
3 EeV the spectrum hardens to s ~ 2.6, a feature known as the ankle. Abo ve 30 EeV one 
observ es a flux suppressio n that has been iden tified with the GZK cut off j Abraham et alj 
l2010ah A recent proposal (lAloisio et al.ll201 11) interprets those features as a proton cut off 
around 10 18 eV and another one around 3 x 10 19 eV associated with iron nuclerl. 



5.2 Implications of the maximum energy 

Cosmic rays at energies below 1 PeV are almost certainly galactic in origin, and those at 
energies above 10 EeV are most likely extragalactic, but considerable uncertainty exists at 
intermediate energies. It is unclear at what energy the local cosmic rays turn from being 
predominantly galactic to being mostly extragalactic. The relevance of this uncertainty for 
modeling the sources of cosmic rays is obvious: if the particles in the energy band above the 
knee at a few PeV, or above the iron knee at ~ 10 17 eV, are extragalactic, then considerable 
finetuning is required in matching the galactic and extragalactic components, because the 
spectrum softens at the knee. On the other hand, if cosmic rays up to a few EeV, i.e. up 
to the ankle in the spectrum, are galactic, then no such finetuning is required, but we need 
to identify the sources of EeV-band cosmic rays with objects present in the Galaxy, e.g. 
supernova remnants (SNR), pulsars, etc. This can be difficult, not because the source in 
question would not accelerate particles to high energies, which in fact we observe happening 
in SNR and pulsars, but because it is questionable that EeV energies can be reached. 



This scenario implicitly postulates a very large proton-to-helium ratio in the source 
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In fact, for typical interst ellar magnetic field values, SNR s hock fronts can hardly ac- 
celerate cosmic rays to a PeV dLagage & Cesarskvlll983alll983bl) . Particle confinement near 
the shock is supported by self-generated magnetic turbulence ahead of and behind the shock. 
Various plasma instabilities driven by cosmic rays can contribute to excite the turbulence to 
high levels, although which dominates remains an active topic of research. In the case of 
SNRs, that which has received most atte ntion so far is the so-called s treaming instabil- 
ity seeded by the cosmic-ray net current dWentzelHl974l : ISkillin3ll975l : lAchterberg|ll983l : 



iLucek & Belli 120001; iBell & Lucekll200ll. e.g. ). and more recently its non-resonant coun 
terpart dBellll2004 120051: IPelletier et al.ll2006h . In contrast, relativistic shocks operating at 



the interface between AGN/GRB flows and the surrounding medium reveal a short pre- 
cursor, which restricts the plasma instabilities to small scale mode s (Medvedev & Loeb 
1999l;|Pohl & ScMickeisdl200d: |Pohl et"aTll2002l: iReville et alJbood IPelletier et al.ll2009l ; 



Lemoine & Pelletied201CLl201 lah . as discussed in detail in the following. Clearly, the ampli 



tude of the turbulence sets up th e pitch-angle scattering frequency and thus the acceleration 
rate dMalkov & Diamondll200ll) . In addition, it also sets the scale for the maximum energy, 
to which a remnant may accelerate particles. Although analytical and numerical estimates 
sugg est that cosmic ray s can very efficiently drive magnetic turbulence ahead of the shock 
(e.g. lBvkov etalEoT l|), s o the turbulent magnetic field may be muc h larger than the homo- 



geneous interstellar field dLucek & Belli l2000l : IBell & Lucek1l200lh . large increases in the 
magnetic field streng th do not necessarily tra nslate into a significant increase in the maxi- 
mum particle energy dVladimirov et al.ll2006h . 

For relativistic sources such as AGN or GRB, relativistic shock acceleration can be in- 



voked, but even there certa i n limitations arise (e.g. Gallant & Achterbergj 19991 : lAchterberg et al.l 
l200ll ; |Pelletier et al.ll2009l ; iLemoine & Pelletieril2011bl : lEichler & PohKOllt) . To be shock 



accelerated, a particle that has crossed the shock toward the upstream must be overtaken 
again by the shock. Assuming the shock moves at Lorentz factor fj, the particle must have 
been deflected (by gyration or scattering) through an angle A 8 > 1 / jS^rj while residing 
upstream. This deflection must be accomplished within a time At at least of order Rs/fcc, 
where Rs is the shock radius at which the shock once again overtakes the particle. The factor 
1 /Ps in A 9 arises for subrelativistic shocks on account of the small incremental energy gain 
per shock crossing. The particle must cross the shock ~ 1 //3s times to double its energy, 
with each crossing requiring at least a significant fraction of a gyroperiod. 

The fastest possible deflection is provided by undisturbed gyration in magnetic field 
oriented perpendicular to the shock normal, for which the angle between particle momentum 
and shock normal, 6, increases linearly with time, and the deflection rate AQ / At must then 
obey A9 / At = f5c/r g > c/(RsTs). This inequality sets a maximum energy ii m ax to which a 
particle can be accelerated, because the rate of change of angle presumably decreases with 
particle energy. As r g = pc/ZeB = f5E/ZeB, this corresponds to a maximum energy, in the 
limit of relativistic particles, of 

E m . dx =ZeBR s r s (9) 

Such regular deflection only occurs in a magnetic field that is coherent over scales larger 
than the path length of the particle. In the short precursor of a relativistic shock, this restricts 
B to the background, undisturbed magnetic-field value, and it thus limits the maximal en- 
ergy to a rather small value for typical interstellar-medium conditions. To be noted is that 
Equation|9]may provide a rather academic limit, because in a relativistic shock a perpendic- 
ular magnetic field in the downstream region renders acceleration very inefficient. Also note 
that the magnetic-field strength, B, is supposed to be measured in the upstream frame of the 
shock, i.e. it is the ambient field in the source frame in, e.g., an extrenal shock of a GRB. 
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Note that scatter-free gyration cannot in general confine a CR particle to a subrelativistic 
blast wave in all three dimensions. The particle generally drifts off to the side after gaining 
the potential difference fcEBRs in energy. Some scattering is required which will reduce 
£ max . Also note that we have neglected both adiabatic losses and drift to the periphery of 
the shock front and assumed that being overtaken by a spherical blast wave is sufficient for 
further acceleration. 

For relativistic shocks, escape through the lateral boundaries does not provide a stringent 
constraint on the maximal acceleration energy unless sideways expansion of the blast takes 
place: as viewed in the shock front rest frame, the particle is confined if its gyration radius 
r g.0|sh < ^i_> w i m R± the lateral extension of the shock front. Since r g ^ sb — r g.o/r s ~, with 
r g the upstream gyroradius in the background field, confinement leads to E nYdx < r<?R±eB, 
which is not as restrictive as the previous expression if R± > Rs/r$. In that limit, sideways 
expansion of the blast is negligible and the overall dynamics resembles that of a spherical 
blast wave. 

What of often invoked magnetic-field amplification by cosmic-ray induced instabilities? 
The growth of plasma instabilities in the precursors of shocks is inevitable, and therefore in 
a realistic situation we cannot expect undisturbed gyration in perpendicular magnetic-field. 
In fact, random scattering is required if the large-scale magnetic field is oriented parallel 
to the shock, because otherwise particles could not return to the shock. Therefore, random 
scattering in small-scale fields will make acceleration at parallel shocks faster, and thereby 
increase the maximum energy, in particular at nonrelativistic shocks. 

The scattering mean free path can be written as A ~ c/Dqq, where the angular dif- 
fusion coefficient is given by Dgg = 8G 2 /8t ~ (eB rms / jSFmc) 2 1/ jSc where 8t ~ I /fie is 
the scattering coherence time, over which the particle scatters by an angle 86, and I is th e 
coherence length of the magnetic field (e.g. lEichler & Pohll I2O1H ; IPlotnikov et alj|201 lh . 
At sub-relativistic shock waves, one must now impose r g /l < PsRs> with r g = pc/ZeB lms . 
For relativistic shocks, the condition for the particle to suffer a rms deflection 1 /J5 over a 
timescale Rs/c reads r g /l < r^Rs, so that these two equations can be combined into 



Emax 

<ZeB ms (p s lR s y r s (10) 



which is less than the previous expression when / <r g , Such small-scale fields are expected 
in the precursor of relativistic shocks, which cannot exceed r g ,o/F s 3 , although there is then 
an ambiguity related to the reference frame of the small-scale magnetic inhomogeneities; for 
simplicity, we have assumed here that these magnetic inhomogeneities are at rest in the up- 
stream plasma. Thus, provided some large-scale perpendicular magnetic field exists, simply 
tangling the field on small scales, /, does not necessarily raise £ max - Magnetic-field ampli- 
fication enhances the maximum energy only if it increase s Bj. ms l. Note that the expression 
Emax = ZeBRs, often taken from Figure 1 of iHillasI J 19841) . is consistent with equations l[9} 
or d 1 Ob only if jS, and F s are both of order unity. We stress again that the limits described here 
may not be reached at a real shock. Leakage from the precursor and the conditions down- 
stream must also be considered when evaluating the maximum energy and the acceleration 
efficiency. 

Equations [9] and [10] suggest that the sources of UHECRs are likely systems involving 
relatistic shocks. 
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5.3 Sources of UHECRs 

Besides reaching the required particle energy, the sources of UHECR must also be powerful 
enough to provide the source luminosity needed to sustain the local flux of UHECR. Possible 
source candidates of UHECRs are active galactic nuclei (AGNs) i Biermann & Strittmatterl 



19871: lTakaharalll990l : IRachen & BiermannI [l993l: IPe'Er et al.ll2009h . clusters of galaxies 



dKang et alj II 19971) Magnetars dAronsI 12003: iMurase et alj 12009^ 7 and gamma-ray bursts 



(GRBs) dWaxmar] ll995l : IVietri|[l995l : lMurase et alj|2006h . Depending on the model, these 
sources may also dominat e the energy range around 10 18 eV feerezinskv et al.l2006l : IWang et al.l 
l2007l : lMurase et al.ll2008h . 

One can constrain the a cceleration capabilit i es of various sources through the magnetic 
luminosity of these sources dNorman et al .1 19951 : IWaxmanll2005l : |l emoine & Waxman 2009, 
e.g.), as follows. Let us assume that acceleration takes place in an outflow at radius r moving 
with possibly relativistic velocity /3 (and Lorentz factor F) towards the observer, so as to 
benefit from Lorentz boosting. We assume that acceleration proceeds with an acceleration 
timescale f acc = J^r g /c in the comoving frame, with si > 1. Then the maximal energy 
at acceleration is at least bounded by the condition r acc < r/(F/3c), which means that the 
acceleration timescale must be shorter than the comoving age of the outflow. This limit can 
be rewritten in terms of the maximal energy in the observer frame, £ max and in terms of 
the magnetic luminosity of the source, Lg = r 2 @ 2 r 2 f3cB 2 /A as calculated in the source rest 
frame in terms of the jet half opening angle © and comoving magnetic-field strength B: 

< 10 20 eV^r-^-'p-^ZL 1 ^ (11) 

with L B ^s = Lb/10 45 erg/s. One can check that this bound remains robust in the small © 
limit, meaning @F — > for which side escape becomes important, and in the small jS 
limit. This bound indicates that rather extraordinary luminosities are required to accelerate 
particles to ultra-high energies, under rather general conditions, although the bound de- 
pends on the charge of the particle. For instance, if one derives the magnetic luminosity of 
blazars through a leptonic modelling of the spectral energy distributions, one concludes that 
only the rare flat spectrum radio quasars with jet powers > 10 44 . . . 10 46 erg/s can accelerate 
protons to ~ 10 20 eV, while other Bl Lac and TeV blaz ars (FR I analogs) with jet powers 
~ 10 40 . . . 10 44 erg/s appear limited to ~ 10 18 - 10 19 eV dLemoine & Waxmanll2009h . From 
this point of view, more compact sources such as GRBs and magnetars appear favored. For 
instance, a GRB of apparent isotropic luminosity 10 52 erg/s with F ~ 100 may produce par- 
ticles with energy as high as Z x 10 21 eV for a magnetic conversion factor ^b = 0.01. As 
discussed in Sec l6.3l mildly or sub-relativistic shocks in a relativistic flow are more efficient 
accelerators of protons than ultra-relativistic shocks and are excellent candidates for being 
sources of UHECRs, owing to the magnetic-field amplification at shocks. 

The paucity of FR 2 radio-galaxies in the GZK sphere (radius ~ lOOMpc) capable of 
accelerating protons to ultra-high energies might be compensated by the acceleration of 
heavier nuclei in the less powerful and more numerous FR I radio-galaxies. In particular, 
Ptuskin and collaborators have shown that if radio-galaxies inject of light to heavy elements 
with a rigidity dependent maximal energy following Eq. (lilt . Lb being related to the radio 
luminosity, accounting for th e radio-luminosity f unction, one could explain rather satisfac- 
torily the observed spectrum dPtuskinetalfeOllI) . It is also intriguing that the Pierre Auger 
Observatory reports an excess of events in the direction of the nearby radio-galaxy Cen A 
(although this latter happens to lie in front of one of the largest concentrations of matter in 
the GZK sphere, the Centaurus supercluster). However, it would be very difficult to under- 
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stand the observed pattern of anisotropy i f one assumes that the high est energy particles are 
heavier than hydrogen in such scenarios dLemoine & Waxmanll2009f) . 

Besides the actual source phys ics, the evolution of sources, the numb er of accelerators 

within a source I Aloisio et al .l2007h . and the variation of source properties dKachelriefi & Semikozl 
l200d : lBerezhkoil2008h will also shape the local spectrum of UHECRs. 

Many properties of UHECRs can be impacted by their propagation in intergalactic 
space, such as their composition through photo-desintegration or their spectrum through 
cascading via photo-meson and photo-pair production, but it is difficult to disentangle the 
propagation effects from the results of physical processes operating inside the sources of 
these particles. 

Estimating the source luminosity using observed quantities and the known population 
statistics of the sour ces in question is subject to considerable uncertainties. As an example, 
lEichler et al.l d2010h have recently estimated the local UHECR source luminosity, assum- 
ing all particles above the ankle at 4 EeV are extragalactic, and compared that with the 
observed gamma-ray prod uction rate of all GRB. In contrast to earlier studies dWaxmanl 
l2004l : lLe& Derme J l2007h . not the MeV-band gamma-ray fluence was used, which likely 
represents a thermal pool, but the GeV-band emission observed with Fermi-LAT, which 
measures the non-thermal tail of the energy distribution in the GRB primary charged parti- 
cles, which, if hadronic, is the part that could contribute to the UHECR flux. It turns out that 
the UHECR source luminosity is more than a hundred times higher than the total GeV-band 
photon output, which places severe constraints on UHECR models involving GRB . If one 
posits that the Galactic to extra-Galactic transition takes place at ~ 10 19 eV and the MeV 
gamma-ray fluence traces the nonthermal parti cle population, the particle output of GRBs is 
more commensurate with their photon output d WaxmarfeO 1 dh . 

For each source class, one can also estimate a luminosity function, that is the differ- 
ential source density needed to integrate the contribution of the sources over cosmological 
redshift. While only nearby (< 200 Mpc) sources may actually contribute to locally observ- 
able GZK-scale UHECRs, the interaction products of the particles from all more distant 
sources will feed a cascade of energy that is eventually observable as a component of the 
extragalactic gamma-ray and neutrino background, which are two other cosmic messengers 
that are complementary to the charged particl es. The former h as been recently measured 
with unprecedented sensitivity up to 100 GeV dAbdo et aL Il2008h . and thus provides invalu- 
able constraints on, e.g., the so-called dip models, which assume that essentially all parti- 
cle above about 1 EeV are protons. The redshift of the onset of photopair production with 
the CMB would then n aturally lead to an ankle at the energy where it is indeed observed 
dBerezinskv et alj|2006l) . More precisely, it is the cosmic evolution of the source class in 
question that determines how much energy is fed into an electromagnetic cascade and even- 
tually reappears in the GeV-band background radiation, relative to the UHECR energy flux 
at the ankle. 



5.4 The transition from galactic to extragalactic origin 



An open problem in cosmic-ray astrophysics is at what energy we observe the transition 
from a Galactic to an extragalactic origin of particles. The limit on inferred source power 
per unit baryon mass required to sustain Galactic UHECR in the [4-40] EeV range that is 
imposed by the observed anisotropy limits is smaller by ne arly 3 orders of mag nitude than 
what is required for an extragalactic origin, as calculated in dEichler et al.ll201C ), and it cor- 
responds to the power per unit mass of gamma rays from GRB ( Eic hler & Pohl [201 lh . This 
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is not only confirmation of the hypothesis that UHECR beyond the ankle are extragalactic, 
it also suggest that their sources are systems not persistently present in the Galaxy. Any as- 
trophysical source class, that is capable of accelerating particles to very high energies and 
should exist in the Galaxy, may fall short of accounting for the trans-ankle UHECR, but may 
nevertheless significantly contribute to the observed cosmic-ray flux between the knee and 
the ankle. For example, the numerical coincidence fits the hypothesis of a GRB origin for 
the Galactic component of UHECR, without invoking a much larger unseen energy reservoir 
for GRB. 

The interesting question is the rate with which such sources appear in normal galaxies 
such as the Milky Way. In other words, what is the role of intermittency? Generally, GRBs in 
the Galaxy are expected every million years or so, the exact rate depending on the beaming 
fraction and the detailed scaling of long GRB with star formation and metallicity. Therefore, 
only a small number of GRB can contribute to the particle flux at the solar circle, and their 
relative contribution depends on the location and explosion time of the GRB. Variations in 
the local particle flux must be expected, and neither the particle spectrum from an individual 
GRB nor the spectrum calculated for a homogeneous source distribution are good proxies. 

IPohl & Eichled d201ll) have calculated the time-dependent transport of UHECR in the 
Galaxy, assuming it can be described as isotropic diffusion. They find that intermittency be- 
comes serious if the mean free path for scattering exceeds 100 pc, unless the source rate is 
much higher than 1 per Myr. On average, Galactic long GRB need to contribute only about 
10 37 erg/s in accelerated particles to fully account for the observed particle flux at 10 18 eV, 
assuming a Bohmian mean free path at this energy. UHECR from Galactic long GRB can 
meet the observational limits on anisotropy only if the mean free path for scattering is suffi- 
ciently small. Contributing the observed sub-ankle particles (at 10 18 eV) requires Bohmian 
diffusion if the UHECR are as heavy as carbon. A light composition such as protons or 
helium requires sub-Bohmian diffusion, which is a highly unlikely situation for isotropic 
diffusion. 

Much of the UHECR anisotropy arises from the expected location of long GRB in the 
inner Galaxy. Observations of GRB host galaxies suggest that regi ons of low metallicity 
and active star formation may be the preferred sites of long GRB dLevesque et alj[20ld : 
iLevesquellioTlh . which may skew the galactocentric distribution of long GRB toward the 
outer Galaxy. As there is no power problem with Galactic GRB, it may be worthwhile to 
also consider short GRB. They provide supposedly less power as a po pulation, but they may 
have a very extended spatial distribution in the Galaxy iBergeill2010l) . leading to a reduced, 
but on account of intermittency not disappearing anisotropy. 

These conclusions can be applied with little change to the case of an origin of UHECRs 
in SNRs, assuming very efficient magnetic-field amplification ca n increase their abil ity to 
accelerate particles to energies significantly higher than 1 PeV (e.g. IPtuskin et al.l2010h . The 
spatial distribution in the Galaxy of long GRB and SNR can be expected to be similar, and 
therefore the average anisotropy is the same for both long GRB and SNR. If one combines 
such a galactic component with a dip model, so that the galactic/extragalactic transition oc- 
curs below 1 EeV, Bohm diffusion and a mixed composition of the Galactic component may 
still be viable, given the systematic uncertainties in the measurements. It would be highly 
desirable to improve anisotropy measurement between 0.1 EeV and 1 EeV, and likewise 
better constrain the composition. 
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5.5 UHECR summary 

Recent progress in UHECR research has built on data from new large-scale observatories. 
The interpretation of measurements of the composition, anisotropy, and spectrum of par- 
ticles provides constraining links between these observables, that are further strengthened 
by new precision measurement of, e.g., extragalactic gamma-ray background emission. The 
very low anisotropy observed for EeV-scale particles provides a strong limit on the con- 
tribution of Galactic sources, if the composition is indeed light as suggested by data. The 
anisotropy found above 60 EeV would be difficult to understand if the particles were heavy, 
which is suggested by Auger data, but not HiRes. If these particles were light, dip models 
might be favorable which, however, must be carefully constructed to not overproduce the 
50-GeV-scale gamma-ray background. 

The main obstacle to further progress clearly is the systematic uncertainty arising from 
the interpretation of the evolution of giant airshowers with particle-physics models that are 
extrapolated over at least 1.5 decades in CoM energy from the range teastable with manmade 
accelerator experiments. 



6 Particle acceleration at relativistic shocks 

Strong shocks occurring in astrophysical flows often generate power-law distributions of 
very-high-energy particles. This is the origin of most high-energy phenomena in astro- 
physics. The favored mechanism for the generation of supra-thermal particles is the famous 
Fermi process. It involves with the scattering of high-energy particles off magnetic distur- 
bances that allow them to cross the shock back and forth and thus to gain energy. Many 
studies have been performed in the 80-ties and 90-ties by assuming pre-existing magnetic 
turbulence. However, it turns out that the pre-existing turbulence is generally not strong 
enoug h to account for the accelera tion performa nce. The nontherm al X-ray emission from 
SNRs (C ass am- Chenai et alJl2004h . but see also IPohl et al.l I2005K revealed that the mag- 
netic field is strongly amplified in the vicinity of the forward shock. Recent theoretical stud- 
ies have shown that the penetration of accelerated particles in the shock upstream flow can 
generate magneti c turbulence that reaches a leve l much larger than the intensity of the am- 
bient mean field dBellll2004l : [Pelletier et al.ll2006h . In producing turbulence the cosmic rays 
loose a fraction of its global energy (about 10 percent of the incoming energy) but increases 
the maximum energy of particles (cf. Section l5T2l . The turbulent field can reach an intensity 
of a few hundreds of jiiG, much larger than the value of a few fiG of the ambient magnetic 
field in the Galaxy. 

These results incited similar investigations for relativistic shocks. Very encouraging re- 
sults were obtained around the turn of the century which extended the theory of Fermi pro- 
cess to the case of relativistic shocks and predicted the formation of a power-law energy 
spectrum with an index s = 2.2 — 2.3 and an acceleration time as fast as the Larmor time 



(Bednarz & Ostrowskil 19981 Gallant & Achterberalll999; 


Kirk et al.l2000l:lAchterbera et al l 


12001: Ellison & Double 120021: Lemoine & Pelletiei 2003 


). But disappointment came once 



the effect of the ambient magnetic field had been taken into account, bec ause it inhibits 
the Fermi process even wh en one considers a strong Kolmogorov turbulence dNiemiec et alj 
l2006l ; lLemoine et al.ll2006h . 



In the following, it will be shown how the paradigm of the three interdependent aspects 
of collisionless-shock physics successfully works in the absence of any mean field: structure 
with a partial reflection on a barrier, supra-thermal-particle generation, magnetic-turbulence 
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generation. Then the scattering issue in the presence of a mean magnetic field will be ad- 
dressed and the requirement for circumventing the inhibition effect will be stated. Then an 
unusual fact in astrophysics will be emphasized, namely the necessity of considering some 
unavoidable micro-physics, that turns out to be crucial not only for the relativistic shock for- 
mation but also for making the Fermi process operative and producing high energy particles. 



6. 1 Successful Fermi process at very low magnetization 

The most favored process for the generation of supra-thermal power law distributions is 
the Fermi process at shocks. Under astrophysical conditions the plasma flow that experi- 
ences a shock is supposed to carry a frozen-in turbulent magnetic field which allows particle 
scattering, and thereby permits particles to gain energy at each Fermi cycle, i.e. a cycle 
upstream-downstream-upstream or downstream-upstream-downstream. 

At a non-relativistic shock of speed f3 s = V s /c <C 1, the average gain per cycle is small, 
G = 1 + | '-y- ji s (where r is the compression ratio, that reaches the value 4 when the shock 
is adiabatic and strong). However this is compensated by a large number of shock cross- 
ings; indeed the escape probability (i.e. the probability for a particle to be entrained by the 
downstream flow and to not come back to the shock front) is low, P esc = 4jS s /r; the return 
probability P let is thus large. A power-law distribution of energy is set up with an index that 
is a simple function of the compression ratio, in the non-relativistic case: 

s=l-—^~l + -. (12) 

In G r — 1 

Strong adiabatic shocks provide a particle spectrum with an universal index, s ~ 2, which 
is modified by losses, radiation losses for the electrons, expansion or escape for protons. 
Subsequent to escape, the spectrum is then steepened by the effect of diffusive propagation 
and escape of particles from the Galaxy. 

A sizable fraction of the incoming energy flux is converted into cosmic ray pressure: 

Pet = $ctP«V? with^ cr -0.1 . (13) 

The successive Fermi cycles produce a precursor of supra-thermal particles (mostly pro- 
tons) of large extension (the diffusion length increases with the particle energy) and this 
penetration in the upstream medium (the ambient medium for an external shock) triggers 
MHD turbulence through two types of strea ming instability, one is res onant and has been 
considered for many years (see for instance iMcKenzie & Voeiklll982h. the oth er is non- 
resonant and has been considered more recently dBellll2004fc IPelletier et al.ll2006l) . as briefly 
discussed earlier. That latter case is quite interesting, first because it is a simple and robust 
mechanism based on the supplementary Lorentz force associated with the plasma current 
that compensates the cosmic-ray current, second because it leads to a turbulent field of large 
intensity; indeed this latter can become much larger than the ambient magnetic field. The 
theory indicates that the fraction of incoming energy flux converted into magnetic energy 
can reach <§g ~ jS s , which is a few percent in SNRs, where one defines 

^=| B p^ 2 . (14) 



A very important remark is that the efficiency of the Fermi process depends on the efficiency 
of the scattering process. By the way, the mechanism of Fermi acceleration is a simple 
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Table 2 Comparison non-relativistic shocks and relativistic shocks. 



At non-relativistic shocks 


At relativistic shocks 


weak escape probability 


significant escape probability 


many cycles of weak energy gain 


few cycles of large energy gain 


power law distribution £~ s with s ~ 2 


power law distribution £ _s with s ~ 2.3 


upstream distribution weakly anisotropic 


upstream distribution strongly anisotropic 


partial reflection at shock front 


partial reflection at shock front 


generation of MHD turbulence upstream 


generation of e.m. micro-turbulence upstream 


acceleration time t acc ~ T s /j8/ 


acceleration time t acc ~ t s 



process, but the scattering, that controls the efficiency of the acceleration process, is the 
main issue. 

As for relativistic shocks, there are similarities and some differences with the non- 
relativistic ones, as summarized in Table 1. There are strong arguments that there is a signifi- 
cant generation of magnetic turbulence at the external shock of a GRB jLi & Waxmanl2006h 
and there is an obvious power-law distribution of ultra-relativistic electrons that synchrotron 
radiate, with an index compatible with the theory of the Fermi process at ultra-relativistic 
shocks (s = 2.2 — 2.3). The ambient magnetic field is very low and at firs t approximation can 
be neglected. A remarkable work was published by ISpitkovskvl ( 120081) that fully validates 
the paradigm, combining three fundamental processes: the formation of a collisionless rel- 
ativistic shock front with reflected particles, the generation of magnetic turbulence and the 
generation of a power-law distribution through the Fermi process. This is a PIC (Particles In 
Cell) simulation of the development of a collisionless shock in a pair plasma (electrons and 
positrons) that runs with a Lorentz factor F, of a few tens (F = (1 — jS 2 ) -1 / 2 ). The flow of 
reflected particles interacts with the flow of passing particles leading to streaming- type insta- 
bilities, and the Weibel branch of instability describes the formation of intense small-scale 
magnetic filaments. The relevant scale of the physics is the inertial length (or skin depth) 
S = . The spatial growth of the magnetic micro-turbulence produces a partial reflection 
of the incoming particles, which allows the formation of a shock front, and self-consistently, 
the reflected particles generate the required level of micro-turbulence. Similarly as the non- 
relativistic case, conversion parameters % C r,%B can be defined in the ultra-relativistic case: 

Per = $ a p u r s 2 C 2 (15) 

% = $B Pu r s 2 c 2 (16) 

And the simulations indicate that <§ cl - ~ 0. 1 and E,b ~ 1 — 10%, similarly to the non-relativistic 
case. The supra-thermal spectrum obtained in the simulation is close to the theoretical pre- 
diction with an index s ~ 2.4. Similar results were obta ined later with PIC simula tion in- 
volving a plasma of electrons and ions of (10... 100) m e |S ironi & Spitkovskvll2009l) . 



6.2 Opening phase space with finite magnetization 

Many astrophysical shocks form in a plasma having a significant magnetization. The physics 
becomes more complex with a finite ambient mean field; it is controlled by the important 
"magnetization" parameter a: 

Bj f Blsm 2 6 B 
4izp u r s 2 c 2 Anp u c 2 
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where Bq is the field measured in the upstream flow frame (generally the ambient field), and 
B t f is the transverse component of the mean field measured in the shock frame. Like in 
non-relativistic shocks, the angle of the field lines with respect to the shock normal is very 
important. But whereas most non-relativistic shocks are in the so-called "sub-luminal" con- 
figuration, i.e. that the angle Gg is not too close to 90° and thus particles can flow along the 
field lines, in ultra-relativistic shocks, it suffices that the field angle 6b be larger than l/T s 
to prevent the return of particles to the upstream region. A generic ultra-relativistic shock is 
thus "supra-luminal", and the magnetic field in the front frame can be considered as almost 
perpendicular, because its transverse component is amplified by a factor F,. This field ori- 
entation is a serious hindrance for the development of Fermi cycles. Neglecting for the time 
being any scattering process in a putative turbulence superimposed on the background field, 
the particle kinematics can be described as follows. A particle that enters the downstream 
flow of speed c/3 is dragged by the frozen in magnetic field a nd cannot easily com e back 
upstream; it can be shown that it can come back just one time dLemoine et al ] |2006h . Once 



upstream, it eventually comes back downstream, but in a subset of phase space that does not 
allow it to make a second cycle. Now, it might be thought that a strong turbulence could pro- 
vide efficient scattering allowing it to make several cycles. However, the typical interstellar 
turbulent field with a large-scale coherence length behaves like an ordered magnetic field 
for such particles, because their penetration length upstream (£ p = m p c 2 /r s eBo, measured 
in the co-moving ups tream frame) is much shorter than the coherence length of turbulence 
dLemoine et alj200a ). In self-generated small scale turbulence, scattering might be efficient 
enough to trigger Fermi acceleration, see below. 

The coherence length i c is formally defined as the range of the field correlation using 
the self-correlation function, C(r). For an isotropic turbulent state we can write (it can easily 
be properly modified in the case of anisotropic turbulence): 



: / C(r)dr; (18) 
Jo 



which can be expressed as an integral over the turbulence spectrum, and one finds that for 
a spectrum proportional to k~P , the correlation length corresponds to large wavelengths for 
1 < ft < 2, as is the case of a Kolmogorov spectrum; for < jS < 1, the coherence length is 
in the shortest-wavelengths part of the spectrum. 

Moreover, the expected duration of the cycle would be much shorter than the eddy turn- 
over time of large-scale vortic es. The requirement s for efficient scattering off magnetic tur- 
bulence are quite challenging dPelletier et al.l2009l) . for not only the intensity of the turbulent 
field must be much larger than the mean field, but also the coherence length must be shorter 
than a Larmor radius. When a scattering process develops, phase space is opened for op- 
erating a Fermi process if the scattering frequency is larger than the Larmor pulsation in 
the mean field. Short-scale turbulence leads to a scattering frequency v, °= e~, whereas the 
Larmor pulsation a>i °c e; thus the range of particle energies for which the phase space is 
unlocked and Fermi process operative, is such that e < e JC0 „ = Ze(B 2 /Bq)£ c . 

At high magnetization (say o > 0.03) the shock is formed by generation of an intense 
coherent wave thr ough a Synchrotron Maser Instability due to a resonance with the loop of 
reflected particles dHoshino & Aronslll99ll ; iGallant et al .111 9921) . The electromagnetic wave 
propagating downstream is damped by synchrotron resonance and produces a thermal dis- 
tribution. The wave propagating upstream carries away a fraction ~ 0.1 CT of the incom- 
ing energy in the case of an e + — e~ -plasma; in a p + — e~ -plasma, an electrostatic wake 
field is generated that heats the electrons up to equipartition while slowing down protons 
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dHoshino & Aronslll99lL ls ironi &Spitkovskvll201lh . The formation of a power-law distri- 
bution jHoshindbOoi) " has not been confirmed as far as we know. The final word has not 
bee n given on these issues of course, because the simulat ions have been conducted s o far in 
ID dHoshino & Aronslll99lMGallant et al.lll992h or 2D dSironi & Spitkovskvll2oTTh over a 
limited amount of time. What happens in a more realistic 3D simulation, the dimensionality 

of which should allow more efficient cross-field transport, or on longer timescales, remain 

I II 1 

to be seen (see |Jones et al.||1998l). 

\.cc< Hili:'.: to lSironi & Spifkovskvl d20 1 lh . at lower magnetization, nothing happens ex- 
cept the thermalization of protons (T p ~ 0.2r s m p c 2 ), until the magnetization reaches a very 
low critical value at which the Fermi process starts. As the magnetization decreases, indeed 
the precursor length scale increases, to the point where plasma microinstabilities triggered 
by the suprathermal particle population self-generate a small scale turbulence that can sus- 
tain the Fermi process (Lemoine & Pelletier 2010). Actually, one needs a very low magnetic 
field to obtain an upstream penetration length of supra-thermal particles large enough for 
having a significant interaction of those particles with the incoming plasma and having a 
growth of micro-instabilities. The Fermi process works with the magnetic component of 
micro- turbulence at the inertial scale ~ 8 = c/(O p i. In principle it starts at even smaller 
scale, the inertial scale of electrons, however electrons are efficiently heated by the electric 
component of micro-turbulence and then the precursor becomes composed of electrons and 
protons of similar relativistic mass, like a pair plasma. This is a very interesting outcome that 
simplifies the physics which rapidly evolves towards conditions similar to those occurring 
in a pair plasma. Thus the PIC simulations of pair plasma are also valuable to understand 
the physics of shocks in electron-proton plasmas. Then a distribution function displaying a 
thermal part and a supra-thermal part with a power law is obtained. 

The transition towards the Fermi process is determined by the micro-instabilities that can 
grow when the upstream penetration of reflected particles is long enough. The fast est insta- 
bilities (Buneman instability, Oblique Two-Stream instability, see lBret et al.ll2004h seem to 
essentially pre-heat the incoming electrons almost up to equipartition with protons. How- 
ever, more simulations are necessary to clarify this important point. The generation of mag- 
netic micro-turbulence by the Weibel instability, which is also studied in laboratory exper- 
iments, is thought to be the main ingredient to form collisionless shocks and to produce 
the Fermi process; however this is also under study by PIC simulations. The generation of 
magnetic micro-turbu lence occurs when the mag netization parameter falls below the fol- 
lowing critical value dLemoine & Pelletiedl20T lal) . as confirmed by numerical simulations 
dSironi & SpitkovskvfcOllh : 

a < a cri , = . (19) 

* s 

Numerical simulations show that the level reached by that Weibel turbulence is such that 
E,b = 1 — 10%, which insures shock formation and Fermi process. Then there exists a large 
energy range for particle scattering when a <S <§J. 



6.3 The micro-physics aspect of GRB termination shocks 

The main issue with Fermi processes based on the scattering off micro-turbulence is that 
the scattering frequency decreases as E , The performance of Fermi processes at non- 
relativistic shocks is determined by the scattering off large-scale, say Kolmogorov, turbu- 
lence which is fairly slow (much slower than the Larmor pulsation in the mean field) but 
decreases only as E^ 1 ^. Thus, if we compare the Fermi process at relativistic shocks with 
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the process at non-relativistic shock, this is like the hare and the tortoise: the scattering, and 
thus the acceleration rhythm, at relativistic shocks is very fast at low energy and decreases 
rapidly as energy increases, whereas, at non-relativistic shocks, it is slow at low energy but 
continues at higher energies with a moderate decline of its efficiency. 



6.3.1 Electron acceleration and radiation 

The external shock that drives the afterglow emission of GRBs may give rise to an effi- 
cient acceleration of electrons if the external medium is weakly magnetized, for the reasons 
discussed previously. If electrons thermalize with protons (as reasonably expected), their 
temperature is already very high at the beginning of the afterglow: T e ~ T p ~ Q.2r s m p c , 
which corresponds to a few tens of GeV. Intense short-scale magnetic turbulence develops 
because the interstellar magnetization parameter is very low, a ~ 10~ 9 , whereas the critical 
value CT crit ~ 1CT 6 , with r s ~ 300. 

What kind of radiation can be expected in such small-scale field, much more intense 
than the mean field? This depends on a so-called "wiggler" parameter a: 

g5 fWE N ^5. (20) 
m e c- m e 

This parameter measures the capability of the magnetic force to deflect a relativistic electron 
of Lorentz factor y by an angle 1 /y (which is the reason why y does not appear in the def- 
inition). If a > 1, then the magnetic field produces a single deflection of the electron in the 
emission cone of half angle 1/y, whereas if a < 1 the electron can undergo several wiggles 
in the emission cone. When a is large, the emission behaves like normal synchrotron radia- 
tion in a mea n field, except th at there is no polarization. When a is small, the emission is of 
"jitter" type dMedve dev 2000). Thus the emission caused by shocked and accelerated elec- 
trons at a relativistic shock is "synchrotron-like", and the analysis of the emitted spectrum 
provides a diagnostic of the magnetic turbulence. 

It is quite remarkable that there exists an almost unive rsal energy limit for th e electron 
radiating in the intense small scale field (in agreement with lKirk & Revillell2010l) : 

w «(«^)^ ( J!^)iAW. (2D 



' Oj m e c 1 ' 



l e < e 



The corresponding maximum energy for the photons emitted in the quasi-homogeneous field 
is 

Er „^'*£^~2x<^(JL) J GeV , m 

where CCf is the fine structure constant. The account for magnetic fluctuations of scales 
larger than the synchrotron emission for mation length resu lts in the photon spectra extended 
beyond the limit given by Eq.@2J (see lBvkov e7al]l2012l) . Thus a single synchrotron-like 
spectrum extending up to several GeV, even possibly a few tens, can be expected and is 
in fact compatible with observations. So the performance of relativistic shocks for electron 
acceleration and radiation appears very satisfactory. The conversion factor into radiation is 
<§rad ~ <§b Cr "o r s < Ye >> an ^ at tne beginning of the afterglow <§ ra( j ~ £,b ~ 1 — 10%. 
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6.3.2 Proton acceleration limited by the fast decay of scattering 

Protons are expected to be accelerated at least as efficiently as electrons at ultra-relativistic 
shock waves. However, as mentioned previously, the ultra-relativistic Fermi process appears 
unable to push protons up to energies in excess of ~ 10 17 — 10 18 eV, because the scattering 
time and thus the acceleration time increase with E 2 in the self-generated turbulent field, or 
scale with E, but then in the background unamplified field. 

Using one or the other, one does not find numbers significantly different from the limit 
associated to the mean field discussed in Sec. l5.2l (Eq.[9l>: E mdx = ZTseBoRs — Z x (0.3 • 
10 7 GeV). Thus, although an energy of order 10 16 eV is achieved, which is something, the 
result is far from reaching the UHE-range. 

Precise performances of mildly or sub-relativistic shocks are not yet known and require 
more numerical simulations. However, some reasonable estimates are permitted by extrap- 
olating what we know about the two extremes: non-relativistic and ultra-relativistic shocks. 
The main guess is that we can expect a magnetic-field amplification at shocks with a con- 
version factor <§b = 1 — 10%, occurring in MHD regime without severe limitation due to 
the super-luminal configuration, especially for oblique internal shocks (termination shocks 
in the hot spots of FR2 jets might be super-luminal). These assumptions can be applied to 
internal shocks of AGN jets (in particular in Blazars jets), and to internal shocks of GRBs, 
as already discussed in Sec. 15.31 



6.4 Conclusion and Prospect 

The triangular dependence of collisionless shock structure with a reflecting barrier for a 
part of incoming particles, with generation of supra-thermal particles and the generation of 
magnetic turbulence is a successful paradigm that applies to astrophysical shocks, both non- 
relativistic and relativistic. Numerical and theoretical works are making significant progress 
for both understanding the physics and providing quantitative results useful for astrophysical 
investigations. This includes not only the spectral index and cut off of the distribution of 
accelerated particles, but also the efficiency factors for the conversion into cosmic rays, 
magnetic turbulence and radiation. We have seen only the beginning of this line of study, 
which requires more PIC simulations and new types of hybrid codes involving relativistic 
MHD coupled with PIC codes for cosmic rays. 

The new results that have already been obtained are important. First, the strong am- 
plification of the magnetic field at SNRs received theoretical and numerical support; the 
astrophysical consequences are interesting, especially for our understanding of the Galactic 
contribution of the cosmic-ray spectrum. Secondly, current state-of-the-art PIC simulations 
indicate that the Fermi process does not operate at ultra-relativistic shocks with magnetiza- 
tion of order unity, which is supposed to be a frequent situation in high-energy astrophysics, 
as for instance in FR2 hot spots, in blazars, in pulsar wind nebulae. Such simulations need 
however to be extended both to higher dimensionality and to larger space-time domains 
before a definite conclusion can be reached. In particular, the issue of the non stationarity 
and/or corrugation of the shock front in relativistic regime should be investigated. Also the 
role of magnetic reconnections in the shock vicinity is a very important new topic whose in- 
vestigation is just starting. Whether and how acceleration proceeds in the mildly-relativistic 
regime also remains open for study. Thirdly the radiation processes that operate in most 
high-energy astrophysical sources involve relativistic electrons scattered in an intense short- 
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scale magnetic turbulence; this leads to a renewed interest in the radiation physics with a 
view to use it as a diagnostic of the magnetic turbulent state. 

The new trend in these topics is the important role imputed to micro-physics phenomena, 
which have a direct astrophysical impact. These developments incite interest in several other 
communities, including space-plasma physics, laser-plasma physics, astroparticle, and high- 
energy astrophysics. We live in exciting times. 
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